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ABSTRACT 

We study first structure formation in A-dominatcd universes using large cosmological Y-body/SPH 
simulations. We consider a standard ACDM model and a AWDM model in which the mass of the dark 
matter particles is taken to be mx = 10 keV. The linear power spectrum for the AWDM model has a 
characteristic cut-off at a wavenumber k = 200 Mpc -1 , suppressing the formation of low mass (< 1O 6 M0) 
nonlinear objects early on. The absence of low mass halos in the WDM model makes the formation of 
primordial gas clouds with molecular hydrogen very inefficient at high redshifts. The first star-forming 
gas clouds form at z « 21 in the WDM model, considerably later than in the CDM counterpart, and the 
abundance of these gas clouds differs by an order of magnitude between the two models. We carry out 
radiative transfer calculations by embedding massive Population III stars in the gas clouds. We show 
that the volume fraction of ionized gas rises up close to 100% by z = 18 in the CDM case, whereas that 
of the WDM model remains extremely small at a level of a few percent. Thus the WDM model with 
nix = 10 keV is strongly inconsistent with the observed high optical depth by the WMAP satellite. 

Subject headings: cosmology: theory - early universe - dark matter - stars:formation 



1. INTRODUCTION 

Popular cosmological models based on Cold Dark Mat- 
ter (CDM) predict that the first stars form in low mass 
(~ 10 6 Af©) dark halos at redshifts z rs 20 — 30 when pri- 
mordial gas condenses via cooling by hydrogen molecules 
(Abel, Bryan & Norman 2002; Bromm, Coppi & Larson 
2002). Hierarchical structure formation eventually leads 
to the emergence of a population of early generation stars 
(Yoshida et al. 2003), which may have at least partly 
reionized the Universe soon after the end of the Dark Ages. 
The first year WMAP result of the measurement of CMB 
polarization implies a large optical depth r e — 0.17 ±0.04, 
indicating that reionization could have occurred as early 
as Zrcion ~ 17 (Kogut et al. 2003). The theoretically 
predicted formation epoch of the first stars in CDM mod- 
els thus appears plausible in light of the WMAP result, 
and indicates that an early generation of stars may have 
contributed to reionization. On the other hand, from the 
theoretical point of view, it is intriguing and important 
to ask whether or not such an early reionization epoch is 
compatible in detail with models other than CDM, or even 
with CDM itself. 

Warm Dark Matter (WDM) models predict expo- 
nential damping of the linear matter power spec- 
trum on small length scales (Bardeen et al. 1986). 
The characteristic length scale is given by the free- 
streaming length of the warm dark matter particle as, 
Rf = 0.31(fi x /0.3) 1 / 3 (V0.65) 1 - 3 (m x /keV)~ 115 h^Mpc, 
where mx is the particle mass (Bode, Ostriker & Turok 
2001). Owing to the suppression of power on small scales, 
the abundance of low mass halos in WDM models is con- 
siderably smaller than in CDM models, perhaps resulting 
in better agreement with the observed matter distribu- 
tion on sub-galactic scales (Bode et al. 2001). Various 



constraints have been placed on the mass mx- Based on 
an analysis of the clustering properties of the Ly-a for- 
est, Narayanan et al. (2000) conclude that the lower limit 
on the mass is mx — 0.75 keV. Dalai & Kochanek (2002) 
used lensing statistics and found that models with mx < 5 
keV are incompatible with the observed abundance of sub- 
structure in distant galaxies. Barkana, Haiman & Os- 
triker (2001) concluded that models with mx < 1 keV are 
likely to be ruled out assuming the reionization rcdshift 
Zreion ~ 6. Reionization at an earlier epoch as implied by 
the WMAP data generically requires early structure for- 
mation, and thus may place a more stringent constraint on 
the mass of dark matter mx 3> 1 keV (Somerville, Bullock 
& Livio 2003). 

In this Letter, we explore the formation of the first bary- 
onic objects in a cosmological model in which dark matter 
is warm, rather than cold. Since models with mx <J 1 keV 
appear to be inconsistent with an array of observations, we 
consider a model with mx = 10 keV. Although the moti- 
vation from particle physics for elementary particles with 
such an intermediate mass is somewhat unclear (but see 
Kawasaki, Sugiyama & Yanagida 1997; Bode et al. 2001), 
it is important to examine the effect of suppressing the 
linear power spectrum on early structure formation. We 
specifically study the formation of primordial gas clouds 
and of their host halos with dark masses 10 5 — 10 6 M Q 
in a ACDM and in a AWDM universe, and compare the 
results for the two models. While the currently available 
observations do not directly probe structure in the redshift 
range we consider here, future CMB polarization experi- 
ments will ultimately reveal how and when the Universe 
was reionized (Kaplinghat et al. 2003), and thus will be 
able to distinguish the structure formation histories pre- 
dicted by the two models. 
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Fig. 2. — Mass function of dark matter halos at z = 20. The 
solid line is the Press-Schechter mass function computed for the 
CDM model, which agrees remarkably well with our CDM simula- 
tion. The thin histograms show the abundance of all the subhalos 
found in the simulation box. 

and normalize the fluctuation amplitude by setting 
cr 8 = 0.9. We follow Bode et al. (2001, see their Ap- 
pendix A) to set-up the initial condition for the WDM 
model, assuming the dark matter mass mx — W keV. In 
order to avoid spurious clumping in the initial particle 
set-up (see Gotz & Sommer-Larsen 2002), we use "glass" 
particle distributions. Further simulation details are given 
in Yoshida et al. (2003, hereafter Paper I). Both of the 
simulations employ 2 x 324 3 particles in a cosmological vol- 
ume of 1 Mpc on a side. The mass per gas particle is then 
160 Mq, whereas the mass of the dark matter simulation 
particles is 1040 Mq. In Paper I, we carried out numer- 
ical convergence tests using higher resolution simulations 
and concluded that the mass resolution adopted here is 
sufficient to follow the cooling and collapse of primordial 
gas within low mass (~ 10 6 Mq) halos. 
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Fig. 1. — The projected gas distribution in the CDM (top) and 
the WDM (middle) simulations at z = 20. Figures in the bottom 
portion show the distribution of dark halos with mass greater than 
1O 5 M for the CDM (left) and for the WDM (right) model. 

2. THE 7V-BODY/SPH SIMULATIONS 

We use the parallel Tree-PM/SPH solver GADGET2 in 
its fully conservative entropy form (Springel & Hernquist 
2002). We follow the non-equilibrium reactions of nine 
chemical species (e~, H, H + , He, He + , He ++ , H2, H^", 
H~ ) using the reaction coefficients compiled by Abel et al. 
(1997). We use the cooling rate of Galli & Palla (1998) 
for molecular hydrogen. We study both CDM and WDM 
models with matter density Qq — 0.3, baryon density 
Sib = 0.04, cosmological constant — 0.7 and expan- 
sion rate at the present time Ho — 70km s _1 Mpc _1 . We 
set the index of the primordial power spectrum n s = 1 



3. RESULTS 

Figure 1 shows the projected gas distribution at z = 20 
for the two models. Panels in the bottom portion show 
the distribution of dark matter halos in each simulation. 
The effect of the exponential cut-off in the initial power 
spectrum is evident in Figure 1; the gas distribution in the 
WDM model is much smoother than in the CDM case. 
Also, the abundance of low mass halos crucial for the for- 
mation of primordial gas clouds is significantly reduced in 
the WDM model. We locate the dark matter halos by first 
running a friends of friends (FOF) groupfinder with link- 
ing parameter b = 0.164. We discard halos which consist 
of fewer than 100 particles. We then use the SUBFIND 
algorithm (Springel et al. 2001) which identifies gravita- 
tionally self-bound sets of particles that are at a higher 
density than the smooth background. We carry out the 
latter step, because, particularly in WDM models, fila- 
mentary structures tend to be identified as halos and such 
objects often contain many gravitationally bound "subha- 
los" (Knebe et al. 2002). Using this two-step method, we 
can robustly identify gravitationally bound objects in our 
simulations. We compare the mass function of the dark 
halos in Figure 2. There we also show the abundance of 
all the subhalos identified in the simulation box by thin 
histograms. The mass function for the CDM simulation 
is well-fitted by the Press-Schechter mass function (solid 
line), as found by Jang-Condell & Hernquist (2001). We 
compare this analytic mass function for the CDM model 
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Fig. 3. — The number of star-forming gas clouds as a function of 
redshift. 

with the result of the WDM simulation (right panel) . The 
difference is nearly two orders of magnitude at a mass 
scale of 1O 5 M . Note, however, that the difference in the 
mass function between the two models becomes smaller 
on larger mass scales, confirming that the suppression of 
the linear power spectrum in the WDM model affects only 
small mass scales. 

In Figure 3, we plot the number of gas clouds against 
redshift. We define groups of cold (T < 500K), dense 
(riH > 500cm -3 ) gas particles as "gas clouds". In order to 
locate dense gas clumps, we again run a FOF groupfinder 
with b = 0.05 to gas particles. From each group we discard 
gas particles which do not satisfy the above criteria. We 
then identify the group as a star-forming gas cloud if the 
cold gas mass exceeds Mj cans = 3OOOM . The first gas 
cloud is identified in this manner at z — 28 in the CDM 
model, whereas it is much later at z = 21 in the WDM 
model. The total number of gas clouds in the simulated 
volume differs by about an order of magnitude in the red- 
shift range plotted. At z = 20, we identified 26 gas clouds 
in the CDM model, and there are only 2 gas clouds found 
in the WDM case. The corresponding numbers at z = 17 
are 66 and 4 for the CDM and WDM models, respectively. 
Note that the number of gas clouds does not represent the 
true abundance of the first stars, because our simulations 
do not include all feedback processes from the first stars. 

An important question is how reionization by an early 
generation of stars proceeds at high redshift in the two 
models. To address this question, we carry out radiative 
transfer simulations for a specific model of star formation, 
as follows. The first stars formed out of a chemically pris- 
tine gas are likely to be very massive (Abel et al. 2002; 
Omukai & Palla 2003). As in Paper I, we assume that 
each gas cloud forms a single massive star, following the 
usual "one star per halo" assumption for "mini-halos." As 
our fiducial model, we set the mass of a Population III 
star to be 300 M with a lifetime of 3 million years, and 
also assume a constant photon escape fraction / csc = 1. 
Although the high escape fraction might seem unrealistic, 
it may indeed be plausible because the gas in mini-halos 
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Fig. 4. — The volume fraction of ionized regions. 

can be wholly ionized and photo-evapolated by a single 
massive star assuming the gas distribution is reasonably 
smooth (Oh et al. 2001). We then run multi-source radia- 
tive transfer simulations using the technique of Sokasian 
et al. (2003). Briefly, the code utilizes a post-procesed 
gas density field defined on a 200 3 grid by casting multiple 
rays from sources in an adaptive fashion. Photon ab- 
sorption and recombination are computed along the rays, 
and each cell carries its own properties such as clumping 
factor and ionization fraction. The gas evolution between 
two adjacent outputs, most importantly recombination, 
is computed on a cell-by-cell basis using these quantities. 
Note that our ray-tracing simulations employ a one-step 
scheme (Sokasian et al. 2001) in which the gas density evo- 
lution due to radiative feedback is not taken into account. 
In order to mimic strong radiative feedback within HII 
regions, we implemented a "volume exclusion effect" by 
disabling sources if they lie within already ionized regions. 
In practice, we turn on a source only if the ionization 
fraction of its surrounding gas is below 0.05. We compute 
the total volume filling factor of the ionized medium from 
the output of the ray-tracing simulation. Figure 4 shows 
the evolution of the volume filling factor computed in this 
manner for the two models. Initially, there is only a single 
HII region and so the filling factor is very small. As more 
stars are formed, the filling factor rapidly increases close 
to 1.0 by z = 18 in the CDM model, causing complete 
reionization. On the other hand, due to the small number 
of sources, the ionized volume fraction in the WDM model 
remains extremely small, only up to 0.03 at z = 18. 

4. SUMMARY AND DISCUSSION 

The suppression of small scale power in the WDM 
model has a significant impact on the formation of primor- 
dial gas clouds. Hierarchical growth of halos with mass 
10 5 — 1O 6 M is not seen in the WDM model, and gas 
cloud formation is nearly completely suppressed until ha- 
los with mass ~ 1O 7 M collapse at z ~ 20. The global 
star- formation rate thus remains very small at z > 17 (see 
Figure 3), regardless of the details of star formation. Our 
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radiative transfer calculations show that reionization by 
early Population III stars is a very slow and inefficient 
process in the WDM model. To be compatible with the 
observed high optical depth of WMAP, the ionization frac- 
tion in the WDM model must increase rapidly at z < 18. 
As Sokasian et al. (2003) argue, the optical depth can 
be as high as r ~ 0.15 only if a large number of ioniz- 
ing photons are produced in (proto-)galaxies at z < 18. 
Clearly the WDM model we consider here is disfavored, if 
not ruled out, in light of the WMAP results. WDM models 
with mx < 1 keV are likely to be ruled out as argued by 
Barkana et al. (2001) and Somerville et al (2003), whereas 
models with mx = 100 keV would be essentially indistin- 
guishable from CDM models in the context of structure 
formation. 

We now turn to the question of whether or not the 
CDM model is compatible with the observed high optical 
depth. Using the results of our ray-tracing calculation 
combined with that of Sokasian et al. (2003) for Popula- 
tion II sources, we compute r e as a function of redshift. 
We combine the contributions from the two modes of star- 
formation in different redshift intervals assuming that the 
onset of Population II begins exclusively at z < 18.5 and 
that there are no Population III stars since then. While 
this is clearly an over-simplification, it provides a conser- 
vative estimate for the total optical depth. To this end we 
re-run the / e sc = 0.20 model from Sokasian et al. (2003) 
with the initial condition that the IGM was fully reionized 
by z ~ 18.5 and was uniformly heated to a temperature of 
1.5 x 10 4 K. Figure 5 shows r e for the CDM model com- 
puted in this manner. Note the slope of the curve decreases 
at 15 < z < 18, reflecting the decline in the ionization frac- 
tion owing to recombinations at a time when the emissivity 
from Population II sources is still low. Population III stars 
in the CDM model give a At ~0.05 whereas the ordinary 
stellar populations contribute At ~0.09 at z < 15, giving 
a total of t ~ 0.14, in good agreement with the WMAP 
result. In Figure 5, we also show r e for the case with only 
the contribution from Population II sources. The reion- 
ization history for the WDM model would be close to this 
case, with the contribution from Population III stars in 
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Fig. 5. — Thomson optical depth as a function of redshift. We 
compute T e using the combined result of the present work for Pop 
III (z > 18) and Sokasian et al. (2003) for Pop II (z < 15) . 

mini-halos being negligible. In a forthcoming paper, we 
will study extensively a number of models using a more 
detailed prescription of early star-formation. The reioniza- 
tion history could be complex as in the double reionization 
scenario (Cen 2002; Wyithe & Loeb 2003), if a dramatic 
transition between the two modes of star-formation, Pop 
III to Pop II, occurs at 6 < z < 18. While it is yet too early 
to draw a definite conclusion, given the uncertainty in the 
WMAP measurement of the optical depth, future CMB 
polarization measurements such as Planck will probe the 
ionized hydrogen fraction at high redshift as well as the 
total optical depth accurately (Kaplinghat et al. 2003), 
and thus will be able to place a strong constraint on the 
structure formation scenario and on the mass of dark mat- 
ter. 

NY thanks Paul Bode for helpful comments on the initial 
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in part by NSF grants ACI 96-19019, AST 98-02568, AST 
99-00877, and AST 00-71019. 
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